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We demonstrate single-photon counting at 1550 nm with titanium-nitride (TiN) microwave kinetic 
inductance detectors. Energy resolution of 0.4 eV and arrival-time resolution of 1.2 microseconds 
are achieved. 0-, 1-, 2-photon events are resolved and shown to follow Poisson statistics. We find 
that the temperature-dependent frequency shift deviates from the Mattis-Bardeen theory, and the 
dissipation response shows a shorter decay time than the frequency response at low temperatures. 
We suggest that the observed anomalous electrodynamics may be related to quasiparticle traps or 
subgap states in the disordered TiN films. Finally, the electron density-of-states is derived from the 
pulse response. 



Fast and efficient photon-counting detectors at near- 
infrared wavelengths are in high demand for ad- 
vanced quantum-optics apphcations such as quantum key 
distribution [l| and linear-optics quantum computing_2:j. 
Superconducting detectors, including superconduct- 
ing nanowire detectors fsj and transition-edge sensors 
(TES)[4J, show great promise in these applications. For 
example, TES made from tungsten films have shown 
over 95 % of quantum efficiency and photon-number 
resolving power at 1550 nm^Sj. On the other hand, 
microwave kinetic inductance detectors (MKIDs) have 
quickly developed into another major superconducting 
detector technology for astronomical instruments from 
submillimeter to x-ray Q- The main advantages of 
MKIDs are that they are simple to fabricate and easy to 
multiplex into a large detector array. Recently, the AR- 
CRON camcrafs^, a MKID array made from titanium ni- 
tride (TiN) films developed for UV/Optical/NIR imaging 
and spectroscopy 0, has been successfully demonstrated 
at the Palomar 200-inch telescope. These TiN MKIDs 
have already shown good photon-counting and energy- 
resolving capability, but so far they have been consid- 
ered only for astronomy applications. In this letter, we 
describe a single-photon-counting experiment at 1550 nm 
with TiN MKIDs and discuss their promise for applica- 
tion in quantum optics. Another motivation for the work 
in this letter is to use these detectors to study the elec- 
trodynamics and microwave properties of TiN, which is 
a relatively new material for MKIDs. Anomalous elec- 
trodynamics of TiN is discussed and the density of states 
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Nq = 3.9 X 10^" eVVm-^ is derived. 

MKIDs are thin-film, high-Q superconducting micro- 
resonators whose resonance frequency fr and internal 
quality factor Qi (or internal dissipation) change when 
incoming radiation with photon energy above twice the 
gap energy {hv > 2A) breaks Cooper pairs in the 
superconductor 3. The measurements of frequency shift 
and internal dissipation signals are referred to as fre- 
quency readout and dissipation readout, respectively. 
The principle of operation of the MKID, as well as the 
readout schemes, were explained in detail in a recent re- 
view paper [Tol|. 

A recent breakthrough in MKID development is the 
application of titanium nitride (TiN), a new material for 
superconducting resonators that shows a number of ideal 
properties for MKIDs [ll|. TiN films have large normal 
resistivity, high kinetic inductance and low loss. Res- 
onators made from stoichiometric TiN films (Tc ~ 4.5 
K) fabricated at NIST and JPL both show Qi > 10^^ 
|l3|. In addition, Tc of sub-stoichiometric TiN films can 
be adjusted between and 4.5 K by varying the nitrogen 
concentration during deposition. This tuning allows the 
gap energy and recombination time to be engineered for 
a specific application. A lower Tc is desired for photon- 
counting MKIDs, because the detector responsivity in 
terms of the frequency shift per incident photon scales as 

The device used in this work is made from a 20 nm TiN 
film with Tc ^ 1 K deposited on a high-resistivity Si wafer 
(> 10 kft ■ cm). The detector design is adapted from the 
optical lumped-element (OLE) MKID design used in the 
ARGONS arrayQ and is illustrated in the inset of Fig. [T] 
The inductive part of the resonator, where photons are 
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FIG. 1: (Color online) Spectra of the signal (solid lines), noise 
(dashed lines) and the signal-to-noise ratio (dotted lines) in 
the frequency quadrature (blue) and dissipation quadrature 
(green) measured at 100 mK. The TLS noise (red) is inferred 
from the difference of the noise in the two quadratures. The 
resonator shows fr = 5.8 GIfz and Q ~ 13, 000 (internal 
Qi ~ 100, 000). We use a readout power Pfeediine ~ —95 dBm 
that is 1 dB below resonator saturation (bifurcation). Inset: 
the optical lumped-element (OLE) MKID design 



absorbed and detected, has been specially tapered to give 
a uniform current distribution and position-independent 
response. Resonators with design frequency between 4 
and 7 GHz are placed directly next to the center strip of a 
coplanar waveguide (CPW) feedline (without the ground 
strip in between) to achieve coupling Qc ^ 10, 000. 

The detector wafer was initially etched by use of Cl- 
chemistry and showed unexpected high level of frequency 
noise. It was then re-etched briefly for 3 seconds in SFg- 
chemistry to clean off the surface residue left by Cl-etch 
that was recently found to significantly increase the loss 
and noise in TiN resonators [14i|. The detectors are cooled 
in an adiabatic demagnetization refrigerator (ADR) to 
100 mK and measured in the standard MKID readout 
(homcdyne) scheme with a high electron mobility tran- 
sistor (HEMT) amplifier (noise temperature r„ ~ 4 K). 
A single-mode fiber going into the device box is used to 
illuminate the entire chip. Outside the ADR, we use a 
1550 nm laser diode driven by a function generator to 
generate optical pulses with width of 200 ns at a repe- 
tition rate of 1 kHz. The pulse response signal of the 
MKID is projected into the frequency and dissipation 
quadratures, digitized at a sampling rate of 2.5 Ms/s, 
and processed using a standard Weiner optimal filter, to 
produce photon-counting statistics. 

The detector responses in the frequency quadrature 
and in the dissipation quadrature measured at 100 mK 
and averaged over 10000 optical pulses are shown in 
Fig. ISlJa). The pulse (peak) height is about 5 times 
larger in the frequency quadrature than in the dissipa- 
tion quadrature, and interestingly, we also see two differ- 
ent pulse decay times, 65 fis for the frequency quadrature 
and 7 fis for the dissipation quadrature. This anoma- 
lous quadrature-dependent response time is discussed be- 



low. On the other hand, the measured noise spectra (see 
Fig.[T]) is consistent with the two- level system (TLS) noise 
picture 15!]: there is 1// excess noise below 10 kHz in 
the frequency quadrature and no excess noise above the 
HEMT noise fioor in the dissipation quadrature. The sig- 
nal and noise spectra shown in Fig. [1] are used to create 
the optimal filter. Due to the reduced pulse height and 
shorter response time, the dissipation readout, although 
free of excess TLS noise, yields a SNR lower by a factor 
of ~ 4 than the frequency readout. Therefore, we used 
frequency readout for photon counting experiments. 
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FIG. 2: (Color online) A histogram of the optimally filtered 
pulse height (normalized by the template pulse) using the 
frequency readout. A 4-peak Gaussian fit to the data is shown 
by the red curve. Inset: the area in each Gaussian peak 
normalized by the total area (blue) is compared to a Poisson 
distribution of A = 1.16 (red). 

Fig. [5] shows a histogram of the optimally filtered pulse 
height using frequency readout. The first 3 peaks are 
clearly resolved, and correspond to the events of 0, 1, 
and 2 photons being absorbed in the detector. We fit the 
histogram to a model with a superposition of 4 Gaus- 
sian peaks, indicated by the red curve in Fig. [5] The 
FWHM widths of the first 3 Gaussians (uncertainty is 
large on the 3-photon peak) correspond to the energy 
resolutions of AE = 0.27 eV, 0.40 eV and 0.44 eV, re- 
spectively. The baseline, A£' — 0.27 eV ^ agrees with a 
SNR analysis from the signal and noise spectra in Fig. [1] 
We conclude from the 1-photon peak that our TiN pho- 
ton counting detector has achieved an energy resolution 
of Ai? — 0.4 eV, corresponding to an energy-resolving 
power oi R — E/AE = 2 at 1550 nm, which is consistent 
with the energy resolution of TiN MKID demonstrated at 
UV/Optical wavelengths foj and comparable with the en- 
ergy resolution of a TES detector operating at the same 
wavelength [5|. We have also derived an arrival-time res- 
olution of At = 1.2 /is from the statistics of the arrival 
time of the optimally filtered pulses. 

It is expected from the stochastic nature of the photon 
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detection process that the n-photon events should obey 
Poisson statistics. Indeed, the counts in the n-photon 
peak (proportional to and calculated by the area of each 
Gaussian) normalized by the total counts match a Pois- 
son distribution of A = 1.16, as shown in the inset of 
Fig. [2j In other words, our detector detects an average 
of 1.16 photons per laser pulsing event in this experiment. 

The two different pulse decay times at 100 mK suggest 
interesting electrodynamic of TiN. We carry out further 
investigations by repeating the photon counting experi- 
ment at elevated bath temperatures between 100 mK and 
190 mK. As we raise the bath temperature, we find that 
the pulse decay time decreases in the frequency response 
and increases in the dissipation response. At 190 mK, the 
two quadratures show the same decay times of 16 /xs, as 
shown in Fig.[3Ub). To explain this phenomenon, we pro- 
pose a hypothesis that low-energy non-dissipative quasi- 
particle (QP) states exist in TiN that give only a reactive 
response but no dissipative response. At lower tempera- 
ture (100 mK), the non-dissipative QP states are mostly 
vacant. The high-energy QPs (generated by photons) 
quickly (within 7 /is) fall into these non-dissipative QP 
states and then recombine slowly (within 65 /is), leading 
to an anomalous response with different decay times for 
the two quadratures. At higher temperature (190 mK), 
the non-dissipative QP states are mostly filled up with 
thermal QPs. The high energy QPs directly combine out- 
side these non-dissipative QP states, giving a normal re- 
sponse with the same decay time for the two quadratures. 
One of the physical candidates for the non-dissipative 
QP states is quasiparticle traps. TiN films are known 
to be highly disordered, with localized superconductiv- 
ity and spatial fluctuations of the gap parameter^lG]. In 
particular, our sub-stoichiometric TiN with Tc ^ 1 K 
sits at a point where Tc is most sensitive to the N2 
concentrationfll]. As a result, we find from DC mea- 
surement that Tc is highly non-uniform and varies from 
0.4 K to 1.3 K over the entire wafer. It is possible that 
a locally depressed gap can form a potential well and 
trap a QP|16|. On one hand, because QPs in the traps 
are frozen in motion, they do not absorb microwave and 
can be dissipationless. On the other hand, these QPs 
are removed (excited) from the Cooper-pair condensate 
and still contribute to a change in the penetration depth. 
This is because, for example, the penetration depth is 
related to the density of superconducting electrons by 
Al c>c 1/ ^/nl in the London theory [l7j. 

The hypothesis of non-dissipative QP states in TiN 
is further supported by the temperature-dependent fre- 
quency shift and dissipation data that also show anoma- 
lous feature (Fig. [3l^c)). Compared to Al resonators, the 
TiN resonator shows excess frequency shift seen as an 
initial slope in 5f{T) at low temperatures, which de- 
viates from the prediction of the Mattis-Bardeen's for- 
mula. This feature has been reported previously in 
NbTiN resonators, where a modified MB formula with a 
gap-broadening parameter F (Dynes parameter jlS]) was 
able to fit the data|is|. In our case, F/A = 0.01 gen- 



erates a reasonable fit to the data. Further investiga- 
tion shows that the density of excitation states modi- 
fied to include the Dynes' parameter has actually intro- 
duced sub-gap states near the Fermi level, which play 
an important role. Excitations to these sub-gap states 
at T < A/fc, which was originally forbidden, gives rise 
to the low-temperature slope in 5f(T). Interestingly, in 
this model these sub-gap states also give a much reduced 
dissipation response. For example, the real part of the 
complex conductivity cti(T) calculated from a QP den- 
sity of states with/ without the Dynes' parameter, or with 
a wide range of F (as long as F/A <C 1) shows little dif- 
ference. Therefore, sub-gap states are another physical 
candidate for the hypothesis of non-dissipative QP states, 
and the anomalous temperature and pulse response are 
consistent under this hypothesis. In fact, sub-gap states 
deduced from tunneling experiments have been reported 
for TiN films lii. 
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FIG. 3: (Color online) Anomalous response of TiN MKID. 
(a) Pulse response at 100 mK. (b) Pulse response at 190 mK 
(c) Temperature-dependent frequency response ^^''^^}^^^'' and 
dissipation response ^^^'^^^'^^ measured for 20 nm TiN {Tc — 
1 K, a = 1, 7 = 1) as compared to 200 nm Al (Tc = 1.2 K, a = 
0.276, 7 = 1/3). a is the measured kinetic inductance frac- 
tion [2l| and 7 is a film-dependent index that takes values 
of 1 in the thin-film limit and 1/3 in the extreme-anomalous 
limit [23]. Fits using Mattis-Bardeen formula with F = and 
r = 0.01 are indicated by the dotted lines and solid lines, re- 
spectively. In all the plots, frequency responses are presented 
in blue and dissipation responses are in green. 
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Despite the anomalous response at lower temperatures, 
we use the pulse response at 190 mK to estimate the 
single-spin electron density of states, Nq, of TiN. This is 
one of the most important material parameters for MKID 
design. From the measured average pulse response of 
Sfr/ fr = 5.6 X 10^® and the average photon energy of 
E = Xhv, we have derived 



Nn = 



VES2 



= 3.9 X 10^° cV^Vm"^, 



(1) 



where r] is the photon-to-QP conversion factor, 
commonly-assumed to be 0.6, T^c ~ 70 /im^ is the vol- 
ume of the inductive part of the resonator; and £'2 is the 
Mattis-Bardeen factor This value of Nq is about a 
factor of 4 larger than the theoretical value previously 
reported 

Our experiment is meant to be a quick demonstration 
and evaluation of TiN MKIDs for use in quantum optics 
applications. Many aspects in the experimental config- 
uration and data analysis can be improved. For exam- 
ple, in our current configuration photons are hitting all 
parts of the resonator without restriction. This may ex- 
plain the small positive shift of the 0-photon peak seen 



in Fig. [2] and may also contribute to the degradation of 
energy resolution from the baseline resolution. This can 
be improved by using a focusing lens or a mask to restrict 
the light to illuminate only the inductive part of the res- 
onator. Also, no attempt has been made in the current 
device to improve the optical coupling or the quantum 
efficiency. In addition to using an antirefiection-coating, 
fabricating the long inductive strip of the MKID on top 
of an optical wave-guide all on a Si wafer might be an in- 
teresting scheme to achieve high coupling efficiency. Our 
TiN MKID has achieved an energy resolution compara- 
ble with that of a TES but is much slower jSj. To further 
improve the energy resolution or speed up the detector, 
a lower-noise amplifier (such as the quantum-limited am- 
plifiers currently under development (23i [2^ ) is probably 
required. This is because the detector noise is already 
limited by the HEMT noise in both quadratures at the 
main Fourier component of the pulse 10 kHz), as can 
be seen in Fig. [1] 

We thank T. Klapwijk, T. Norguchi, J. Martinis, G. 
O'Nell and A. Lita for useful discussions. The TiN film 
was deposited in the JPL Microdevice Lab (MDL), and 
the device was fabricated in the NIST cleanroom. 
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